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Abstract 

Wilson loops with lightlike polygonal contours have been conjectured to be equiva- 
lent to MHV scattering amplitudes in M = 4 super Yang-Mills. We compute such 
Wilson loops for special polygonal contours at two loops in perturbation theory. 
Specifically, we concentrate on the remainder function 71, obtained by subtracting 
the known ABDK/BDS ansatz from the Wilson loop. First, we consider a particu- 
lar two-dimensional eight-point kinematics studied at strong coupling by Alday and 
Maldacena. We find numerical evidence that 1Z is the same at weak and at strong 
coupling, up to an overall, coupling-dependent constant. This suggests a universal- 
ity of the remainder function at strong and weak coupling for generic null polygonal 
Wilson loops, and therefore for arbitrary MHV amplitudes in M = 4 super Yang- 
Mills. We analyse the consequences of this statement. We further consider regular 
n-gons, and find that the remainder function is linear in n at large n through numer- 
ical computations performed up to n = 30. This reproduces a general feature of the 
corresponding strong-coupling result. 
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1 Introduction 



Scattering amplitudes in gauge theories, and in particular in Af = 4 super Yang- 
Mills (SYM), exhibit many beautiful properties which have greatly expanded our 
understanding of their mathematical structure. In |T|, a remarkable relation was 
discovered, expressing the planar two-loop four-point amplitude in Af = 4 SYM in 
terms of the one-loop amplitude. The authors of [l] also conjectured a generalisation 
of this relation for two-loop MHV amplitudes with an arbitrary number of external 
legs. The possibility of iterative relations at higher loop order was investigated in 
[2], and confirmed for three-loop, four-point amplitudes. Furthermore, an intriguing 
conjecture was put forward in (2), expressing an appropriately defined finite part 
of the planar all-loop n-point MHV amplitude as the exponential of the one-loop 
amplitude times the cusp anomalous dimension [3]-|5]. An equation for the all- loop 
cusp anomalous dimension was found in j6], and the strong-coupling expansion of its 
solution presented in [7]. 



In a remarkable development, Alday and Maldacena [8] were able to calculate 
the four-point amplitude in Af = 4 super Yang-Mills at strong coupling using the 
AdS/CFT correspondence. Their result confirmed the structure of the ABDK/BDS 
proposal at four points, with a simple replacement of the weak-coupling expression of 
the cusp anomalous dimension with its strong-coupling counterpart. The prescription 
of [8] maps the calculation of scattering amplitudes to that of a Wilson loop whose 
contour is obtained by gluing the particle momenta following the ordering of the 
particles in the corresponding colour-stripped, planar amplitude. This suggested a 
remarkable duality, according to which MHV amplitudes are equal to expectation 
values of lightlike Wilson loops [8 10 . This duality was tested at one loop for four 



particles in [9], and for any number of particles in 10 . Further calculations performed 
in 11 12 confirmed the duality at two loops in the four- and five-point case. The 



iteration of the parity-even part of the five-point amplitude was checked in 13 
This result was extended to the complete amplitude in (ll), where a cancellation of 
contributions arising from parity-odd terms in the logarithm of the amplitude was 
discovered. 



The BDS ansatz was further studied in 15 , where it was found that, for a par- 
ticular kinematic configuration with a large number of scattered particles, a certain 
discrepancy function - now known as the remainder function - should be added to 
the BDS ansatz in order to reproduce the complete amplitude. The same conclusion 
was also reached for the six-point case, based on perturbative Wilson loop calcu- 
lations in fl6l and the studv of Regee limits in 17 . The remainder function was 



[16] and the study of Regge limits in 
then obtained by an explicit two-loop calculation of the parity-even part of the six- 
point MHV amplitude 18 . In 16 19 the corresponding hexagon Wilson loop was 



computed, and the comparison of the amplitude and Wilson loop remainder functions 
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performed in 18,19 showed that these two quantities were identical in all kinematical 
configurations considered. 



One key goal is to characterise and eventually determine this remainder func- 
tion. An important insight of [12] was that the Wilson loop remainder function must 
be invariant under dual conformal symmetry. At four and five points, the lightlike 
condition on the particle momenta do not leave room for non-trivial conformally in- 
variant cross-ratios, and the BDS conjecture remains intact. Only starting from six 
points one can write down such cross-ratios, and dual conformal symmetry implies 
that the remainder function depends on the kinematic variables only through these 
cross-ratios 
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Additional information about this function was found in 18 , 20 using simple 



and multi-collinear limits of amplitudes and Wilson loops. Furthermore, Anastasiou, 
Spence and the present authors developed in 20 a numerical algorithm to evaluate 



Wilson loop remainder functions with an arbitrary number of edges, based on pro- 
grams developed in 



21 24 for evaluating Feynman diagrams. The results of 



make it possible to compare the Wilson loop remainder function to the amplitude 
remainder function, once the integrals of the latter are known to the required order 
in the Laurent expansion in the dimensional regularisation parameter e. We note 
that recently the parity-even part of an arbitrary two-loop n-point MHV amplitude 



has been written down in 25 in terms of integral functions. Particular limits of the 



two-loop Wilson loop remainder function were also recently studied in 26 , 27 



In a very interesting recent development 28 



29 



30 



a new class of Wilson loops was 
In the AdS/CFT set-up of M, 



computed at strong coupling using results of 
the problem of calculating the amplitude at strong coupling is equivalent to that of 
minimising the area of a string worldsheet bounded by a lightlike, polygonal contour 
living on the boundary of the AdSs space. In 29 , Alday and Maldacena considered a 
more tractable situation where the contour of the Wilson loop is embedded in a M. 1,1 
subspace of the boundary of AdSs, with the worldsheet living in an AdS3 subspace 
of AdSs. The simplification achieved in this way is quite dramatic - the number 
of independent cross-ratios for a 2n-gon contour is reduced to 2n — 6, where dual 
conformal symmetry has been used to map the position of three cusps at infinity. 
This eliminates six out of the 2n parameters needed to fix the position of the 2n 
cusps. Hence, at six points there are no free cross-ratios and the remainder function 
is a constant. Starting from n = 4, or eight points, we have non-trivial octagon 
kinematics determined by two cross-ratios. In 29 , the explicit form of the octagon 



remainder as a function of these two cross-ratios was presented in terms of a simple 
one-dimensional integral (up to constant terms). 



The original motivation of this paper was to evaluate at two loops the octagon re- 



mainder function in the kinematics of |29 using the techniques of 20 and to compare 



2 



with the strong coupling expression found in 29 . It is far from obvious that there 
should exist any simple relation between the strong and weak coupling limits of the 
remainder function. However, the results we find suggest that, in fact, the remainder 
functions at weak and strong coupling are linked by a simple linear relation, 

n^( Uij ) = A^n 8 AM ( Utl ) + b® , (i.i) 

where 1Z% (uij) is the two-loop remainder function, and 1Z AM is the strong-coupling 



remainder evaluated in 29 ; the cross-ratios are denoted by u^. Using simple- and 
multi-collinear limits, we will further argue that the constant B^ can be absorbed 
into a redefinition of the strong-coupling remainder function, which in practice leaves 
only one coupling-dependent quantity, A^ 2 \ in the relation linking strong and weak 
coupling. We should stress that our calculation is based on numerical evaluations, 
and is therefore correct only within errors. 

If the correspondence does hold, it suggests the following general matching relation 
between the all-orders perturbative remainder function, lZ(a), and the strong-coupling 
expression, TZ stTOQg , 

K(a) = A(a)TZ stTong , (1.2) 

where a is the 't Hooft coupling, and the function A (a) has a perturbative expansion 
starting at two loops and is independent of n and of the kinematic variables (which 



are all contained in 7^ str0 ng)- If (1-2) is correct, it should be very interesting to find 
the physical meaning of the proportionality constant A(a). Since the BDS part of the 
amplitude already has a universal structure at weak and strong coupling, combined 



with (1.2), this implies the universal structure of the Wilson loop (MHV amplitude) 
itself. 



The second special class of kinematic configuration addressed by 29 is that of 
regular 2n-gons. In this case the remainder function is just a number, and one can 
study it as a function of n. Our main result here is that, at large n, the regular 
remainder function at two loops is linear in n. In this case we do not perform a full 



comparison of our weak-coupling results to the strong-coupling results of 29 , where 



only one part of the complete remainder function was presented. It would clearly 



be of great interest to check whether a relation such as (1.2) still holds. This could 



lend further support to the conjecture that (1.2) is valid not only for the special 
octagon kinematics embedded in K 1 ' 1 , but for all ra-gon Wilson loops, and for the 
corresponding n-point MHV scattering amplitudes in M = 4 SYM. If this is correct, 
the remainder function is universal. 

The rest of the paper is organised as follows. In Section 2, we review some basic 
facts about the Wilson loop remainder function, in particular its precise definition 
and its behaviour in (multi-)collinear limits. Section 3 is devoted to the calculation 
of the octagon remainder function at two loops. We begin by first reviewing the 
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(l + l)-dimensiorial kinematics introduced in 29 , and map it to conformally equiv- 



alent (2+l)-dimensional configurations where all the positions of the cusps are at 



finite distance. Thanks to dual conformal invariance of the remainder function 12 
all these configurations are equivalent for the sake of calculating this quantity, how- 
ever our (2+l)-conflgurations are more convenient for numerical simulations. We 
then compare the results at strong and weak coupling, presenting evidence in favour 
of (1.1), as well as estimates of the parameters and appearing in that rela- 
tion. In Section 4 we analyse the consequences of the universality conjecture for the 
remainder function, and discuss how to reabsorb the constant into a redefinition 
of the strong- coupling remainder function. Section 5 contains our results for regular 
2n-gons. Finally, we present our conclusions in Section 6. 



2 Definition of the remainder function 



We expand the Wilson loop in powers of the 't Hooft coupling a a^j] 



(W[C n ]) :=l + J2 alW n l) ■= expJ2 alw n 

1=1 1=1 



(2.1) 



In particular, 



w 



(i) 



w 



(2) 



1)\2 



(2.2) 



We note that times the tree-level MHV amplitude is equal to the one-loop MHV 
amplitude first calculated in 31 using the unitarity-based approach [32], up to a 
regularisation-dependent factor. Wn^ was calculated in 9 for n = 4 and in 
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for 



arbitrary n, while Wn^ was evaluated analytically for n = 4 in [ll], and numerically 
for n = 5, 6 in |l2|[l6(|l9||20| . Results for n = 7, 8 were also presented in [20] . 



We define the remainder function 1Z^ L for an n-sided Wilson loop via 



w 



(2), 



(2) 



iWL 



where e < 0, and 



f(2) I \ ._ iW i tV-) e i fW 
JWLlW JO i il,WL fc ' J2,WL t ■ 



e(2) , ,(2) 



(2) 2 



(2.3) 
(2.4) 



/q 2 ^ is the same as on the amplitude side, = (2, while /{ 2 WL = G^k = ^C3 
In 20 , the four- and five-edged Wilson loops were cast in the form (2.3) with 



33 



(2.5) 



1 Wc follow the definitions and conventions of 20 , to which we refer the reader for more details. 
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which allowed for a determination of the coefficients wl anc ^ ^wl- The resu hs 
found in (20) , ar(Q 

(2.6) 



(2) / 



and 



As noticed in 



a 



(2) 
WL 



(2.7) 



20 



(2) 

there is an intriguing agreement between the constant C^ L and 



the corresponding value of the same quantity on the amplitude side. 



The statement of the duality between Wilson loops and amplitudes beyond one 
loop can be recast as an equality of the corresponding remainder functions, 



K 



WL 



{21 



A consequence of the precise determination of the constants / 2 W l anc ^ C*wl * s that 



.(2) 



no additional constant term is allowed on the right hand side of (2.8). For the same 



reason, the Wilson loop remainder function must then have the same collinear limits 
as its amplitude counterpart, i.e. 



tjWL . -r>WL 

'Si — f 'Si-i j 



(2.9) 



with no extra constant on the right hand side of (2.9). 



As is well-known [12], the remainder function depends on the kinematics only 
through conformal cross-ratios. In 20 , the following definitions for the cross-ratios 
were adopted, 



U; 



ry 2 ,y. 2 

/ >-» 2 ™ 2 

^ij^i+lj + l 



(2.10) 



where the x« are 't Hooft's region momenta, in terms of which momenta of on-shell 
scattered particles are defined as pi = X{ — Xi+\. The number of these cross-ratios is 
the same as that of two-mass easy box functions, i.e. n(n — 5)/2, if one ignores Gram 
determinant constraints 20 . The inclusion of these constraints reduces this number 
to 3n — 15 [12] . For n = 6 one obtains three cross- ratios in both cases, whereas for 
octagon Wilson loops - the case of principal interest for this paper - one obtains 9 
(12) cross ratios with(out) Gram determinant constraints. In [20] it was also shown 
that a generic cross-ratio, i.e. one of the form x^xf m / (x^x?) can be re-written as an 



appropriate product of the cross-ratios defined in (2.10). 



Finally, we briefly summarise the results of 18 , 20 on the expected behaviour of 



the amplitude and Wilson loop remainder functions under triple-collinear limits. We 



2 The 0(1) and 0(e) coefficients of f^ L (e) had been determined earlier in pd| . 
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consider the case where the three momenta pi, p 2 and p 3 become collinear, 



Pl := X1-X2 = Z\P , P2 = X2-X3 = Z 2 P , p 3 = X 3 - X 4 = Z 3 P , Zi + Z 2 + 23 = 1 

(2-11) 

with P 2 = (px + p 2 + p 3 ) 2 



2 '0. In this limit, one finds 18120 that 



where 



u-i 



1 S12 



1 S23 



1-23S123 l-^lSi 23 (1 - Zl)(l - Z 3 ) 



(2.12) 



(2.13) 



3 Octagon Wilson loops 



3.1 Definition of the kinematics 



We begin by reviewing briefly the special eight-point kinematics introduced in 29 



The starting point of 29 is a "zig-zag" configuration of momenta embedded in 



C 



R 1,3 where three of the eight region momenta Xi which determine the positions of the 
cusps are at infinity, 



Xi = (0,0) , x 2 
1 



X.-i 



1, 



i + r 

x 8 = (0, -00) . 



1 + 

x 5 



x 3 



(3.1) 



1,1) , 



i + x+ i + r 

xq = (00, 1) , xj = (00, —00) , (3.2) 



This configuration depends on two parameters, or more precisely lightcone cross-ratios 
X ± , introduced in [29]. Alternatively one can combine them into a single complex 
parameter m as follows, 



3 27rlm(m) 



X 



-2vrRe(m) 



(3.3) 



matics (3.1), obtaining 



We can straightforwardly compute the 12 cross-ratios (|2.10|) arising from the kine- 

(3.4) 



U15 



X' 



i + r 
1 

i + r 
1 , 



U 2 Q 
U 48 



X' 



1 + 

1 



1 + X" 



1,...,; 
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One of the goals of this paper will be to compare the strong-coupling expression for 
the Wilson loop for the particular kinematics of |29] with perturbative computations. 



Direct use of the (l+l)-dimensional kinematics of (3.1) is awkward for numerical 
evaluations, due to the presence of infinities. In order to circumvent this problem, we 
will employ instead three different (2+l)-dimensional sets of kinematic invariants in 
our perturbative calculations, each of which is conformally equivalent to the set above, 



i.e. gives rise to the same cross-ratios as in (3.4). Since the Wilson loop remainder 



function is dual conformal invariant 12 , in practice one can equivalently determine 



it with any of these three sets. This also provides additional tests for our numerical 
calculations. 

Kinematics A 

Our first choice of (2+l)-dimensional kinematic variables is given by 

Xl = (1/2,1/2,-1), x 2 ' V X 
l + {2 + X -)x + 



x 3 



2 + 2 X +'2 + 2 X + ' 

1 + -(i + x")(i + x + ) 



2(1 + X- + X-X + ) ' 2(1 + X~+ X~X + ) ' 1 + X- + X~X + 
' ' ' x 5 = (1/2,-1/2,-1), 



2 + 2 X -'2(l + X -)' 
x 6 = (-1/2, -1/2, 0) , x 7 = (0, 0, 0) , x 8 = (-1/2, 1/2, 0) 



(3.5) 



It can easily be shown that the corresponding cross-ratios (2.10) are exactly as those 



in (3.4). 



Kinematics B 



An even simpler choice of kinematics, which is conformally equivalent to that in (3.1 ) 



is characterised by taking the cross-ratios of (3.4) and supplementing them by eight 



Lorentz invariants, all taken to be equal to minus one. More precisely, following 20 
we start from the full set of twenty momentum invariants 



X i i+2i X i+4 i+6> X ii+3i X i+4i+li X ii+4 i % 



1,...,4, (3.6) 
and fix the following three- and four-particle momentum invariants as follows: 



x 



i+5 i+8 



-1 



X i i+4 ' 1 



1,...,4 



(3.7) 



The remaining invariants are then determined in such a way that the twelve cross- 



ratios are the same as in (3.4). The computation then proceeds in terms of these 



invariants without reference to specific region momenta. 
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Kinematics C 

The last (2+l)-dimensional kinematics we will use in the following satisfies the addi- 



tional requirement 29 that the spatial projection of the polygonal contour circum- 
scribes the unit circle, see Figure [Un The coordinates of the region momenta are then 



given in terms of two parameters x, y as follows: 

Xi = (tan y, 1, tan?/) , x 2 = (— tan z, cos (z + 2y) sec z, sin (z + 2y) sec z) , 
x 3 = ^tan x, — cos (x + — j sec x, sin ^x + — j sec x 

x 4 = ( — tan ( — ) , — l,tan ( — ) ) , x 5 = (tan ( — ) , —1, — tan 



- tan (J^j , - tan (J^j , x 7 = (tan (^j , tan (J^j , -1 

-tan(^) ,l,-tan(^)) , (3.8) 

where z = 3ir /8 — x — y. 

The cross-ratios computed in this "circular" kinematics coincide with those of 



3 This kinematics is denned in terms of three angles x, y and z — 3tt/8 — x — y. For this polygon 



to circumscribe a circle, the three angles must be positive. More generally, the kinematics in (3.8) 



remains valid for any values of x and y (even if the contour does not circumscribe the unit circle) . 
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(3.4), if x and y are related to \ m the following way, 

, -tanx-\/2 + l 

X = —j=- > (3-9) 

V 2 tan x — tan a; — 1 

-tany-v^+l , in * 

X = —?= • (3.10) 

y2tan?/ — tany — 1 

The above equations can be easily inverted to give explicit analytic expressions for x 
(for y) as a function of x + (x~)- 

Note that in the limit m-^Owe have x * 1> an d hence 

U ii+ 4 = r, 1 = 1,..., 4, (3.11) 
ttii+3 = 1, z = l,...,8. 

We will see in Section [5] that these cross-ratios correspond to the case of a regular 
octagon. For the case of our circular kinematics this is also obvious from Figure [TJ 
since in this limit (x,y) — > (tt/8, 7r/8). 



3.2 Results at strong coupling and limits on the kinematics 



In the following we briefly review the strong-coupling result for the octagon remainder 



function in the special two-dimensional kinematics (3.1) introduced in 29 . This 
remainder function was found to be equal to |29 



Kt u (m) = -ilog(l + l/ x + )log(l + X " 



+ Jl + dt H Shlht l og (1 + e -2.| m |coshA (3 12) 

6 tanh(2t + 2z</>) & V ) i \ ) 
where the complex variable 

m := |m|e^ , (3.13) 

is related to the lightcone cross-ratios x* 1 through ( 3.3[ ) j^j The expression for the 



octagon remainder function in (3.12) is appropriate for the first quadrant of the 
complex m plane (<fr = 0, . . . , vr/2), and is defined in the other quadrants by analytic 
continuation. Furthermore, the remainder function is symmetric under — > —<p and 
— > (f) + 7r/2, hence knowledge of this function in the first quadrant is sufficient. 

For general 2n-point amplitudes (or 2n-gon Wilson loops), the remainder function 
in the notation of [29] can be represented as 



= AA^ S + (A™ h + A^) + AlT° ds - (3.14) 



Recall that these arc related to standard conformal cross-ratios via (13.4 
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The first term, AA™ S := ^BDS-iikc-cven _ a bds^ in the cage of the octagon gives rise 



in (3.12), and finally A 



to the log-log term in (3.12); the second term in (3.14) gives 77r/6 plus the integral 



periods 



0. We find it convenient to separate the contribution 



t sinh 



arising from regular 2n-gons (i.e. when x ± 0), and define Af^ 1 := A i Teg2n 
Summarising, for the octagon we have: 



+ AA 



sinh 
2n ' 



AAg DS 



a sinh 



^yjsinh _i_ ^jextra 
8 8 

^periods 



We recall that A^ riods 



2 

57T 

T' 

71 

-T2 + 

0. 



log(l + 1/ X + ) log(l + X 



dt 



I m I sinht 



tanh(2£ + 2iArg(m)) 



log (1 + 



for regular polygons, while A 



extra 
2n 



(3.15) 
(3.16) 

27r|m|coshi ^ 17) 

(3.18) 

for odd n 



29 



Two limits of these expressions will be important for our discussion, and have an 
interesting geometrical interpretation: m — > and \m\ — > oo. The first limit, m — > 
(or x ± ~^ 1) corresponds to the case of the regular octagon. Note that in this limit 
AAf nh + A| xtra 0, as shown in Q, so that 



m — > 



K 



AM 
reg- 



1 



5n 



2 '< 2+ 4 



(3.19) 



In the second limit, as \m\ goes to infinity, the cross-ratios x either vanish or become 
infinite. For example, in the limit where Re(m) — > +oo and Im(m) — > — oo, one has 
X + — > and x~ ~^ 0- m this case one can show that 



K 



AM, 



m) — > 



5tt 

T 



7T 

12 



7vr 

~6~ 



(3.20) 



Intriguingly, this coincides with twice the value of the regular hexagon remainder 
function, 7^ M 



7vr/12 29 , for which one has (cf. (3.15 )-(3.18 )): 



K 



AM 



AA 



BDS 



+ A 



sinh 
reg— 



+ AAg inh + A 



extra 
6 



7vr 

0+-+0 



0. 



(3.21) 



The fact that TZs 2TZq in this limit has a general physical explanation, also valid at 
weak coupling. In general, the large- |m| limits of 1Z$ correspond to various double- 
soft or soft-collinear limits 29 . In particular, the limit considered above, x ± ~~ ^ 0> 



corresponds to a double-soft limit, i.e. a situation where two consecutive edges are 
becoming small. 



This double-soft limit can be regarded as a special triple-collinear limit, where 
(pi,p 2 ,P3) - > [ziP, Z2P, z 3 P) and furthermore — » (as usual, P := pi + p 2 + P3, 
Z1 + Z2 + Z3 = 1). This implies that the eight-point Wilson loop collapses to a six-point 
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Wilson loop plus the corresponding triple-collinear splitting amplitude. As reviewed 
in Section |2j the latter is nothing but the six-point remainder function 18 , which 
for the special kinematics we are considering is just a constant. Hence, we find that 
for 0, 7^8 — > 7^6 + TZq, where the second 7£ 6 is precisely the contribution of the 

triple-collinear splitting amplitude. In summary we have, 



m — > oo 



211 



AM 
feg— 6 



7vr 

~6~ 



(3.22) 



In order to visualise the strong-coupling remainder function (3.12), and to compare 



it with our result for the two-loop remainder function, we will plot both of them in 
the following section in Figure [2] as functions of Re(m) and Im(m). 



We shall now define a rescaled octagon remainder function at strong coupling, 
which can be directly compared to our weak-coupling results in the next section. To 



this end, we first shift 1Zg M (m) as defined in (|3.12) by its value at m = 0, and then 



divide by the maximal variation 7Z 8 



AM/ 



(0) - K^ioo). Using the limits in <KM and 



(3.22), we get 



nt M (m) 



m) 



-I? AM 
/v reg-8 



K 



AM, 



Hog 2 2 



5tt 

4 



-nAM __ o-pAM 
/v reg-8 ^^reg-e 



\ log 2 2 + 



(3.23) 
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Figure 2: The octagon remainder function in the special kinematics considered in 29 



as a function of Re(m) and Im(m). The figure on the left represents the strong- 
coupling result 1Z$ M derived in ( |3.12[ ) ; while on the right we plot our two-loop re- 
sult fti 2) . 



We now present our results for the octagon remainder function at weak coupling, 
and compare the perturbative remainder function to its strong-coupling counterpart. 
We have computed the Wilson loop at two loops using the general set-up developed 
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in our earlier work 20 , to which we refer the reader for a discussion of the numeri- 
cal routines employed. The remainder function is then extracted from the two-loop 
Wilson loop using (2.3). 



In order to compare our results to the remainder function evaluated at strong 
coupling we are free to use any conformally equivalent kinematic set, as the remainder 
function is dual conformal invariant [12]. As mentioned earlier, for practical purposes 
we will use three different choices of kinematic variables, A, B and C, all of which 
are conformally equivalent to the kinematics of |29], as described in Section 3.1 As 



a preliminary question, one can check dual conformal invariance by comparing the 

(2) 

value of the two-loop remainder function 7Z S at a fixed set of cross-ratios, realised 
through the three different kinematics A, B and CM We have confirmed in a large 
number of cases that this is indeed the case, within the numerical errors. In Table [T] 
below we quote one particularly important point, namely m = 0, which corresponds 
to a regular octagon]^] 



kinematics A 


B 


C 


/v reg-8 


-5.5275 ±0.0008 


-5.52745 ± 0.0008 


-5.5284 ± 0.0007 



Table 1: In this table we list our results for the Wilson loop remainder function 
corresponding to a regular octagon contour. The cross-ratios corresponding to this 
case, quoted in (3.11), are realised using three different kinematics, namely the sets 
A, B, C. 



To begin with, we show on the right hand side of Figure |2l a two-dimensional plot 
of our result for 1Z\ as a function of the two real variables Re(m) and Im(m). For 
comparison we show the strong-coupling remainder function on the left. Already at 
this stage one can note an obvious similarity in the shape of these two functions. In 
the following, we will quantify this similarity in a more precise way. 



To address this point, we begin by making the preliminary observation that the 



geometrical interpretation of the m — )■ and |m| — > oo limits discussed in Section 3.2 



is obviously valid at strong as well as at weak coupling. Thus, (3.19) and (3.22) hold 
also in perturbation theory (at any loop order): 



\m\ — > oo : 
m — t- : 



271 



(2) 

reg— 6 ' 



K 



(2), 



m) 



/v reg-8 ) 



(3.24) 

(3.25) 



and only the specific values for the regular octagon and the hexagon differ between 
5 This provides a useful check of the validity of our numerical routines. 

6 This data point will be used later in order to define the rescaled two-loop remainder function, 



which can be directly compared to the corresponding strong-coupling result in (3.23) 
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the two regimes. At two loops, we get for the regular hexagorQ 



K 



(2) 



-2.707 ±0.002 . 



It was observed in [20] that (3.26) is compatible with 



K {2) -v -— 



(3.26) 



(3.27) 



For the regular octagon we will take as a reference value the following (see Table 0: 

-5.5284 ±0.0007. (3.28) 
In order to compare our weak-coupling result to the strong-coupling calculation of [29] , 



7?( 2 ) ~ 

/v rcg-8 — 



we introduce a shifted and rescaled remainder function, similarly to (3.23) of the 
previous section, 



Tlf(m) :-- 



nf\m) 



K 



(2) 
reg- 



-ip(2) _ 9-77(2) 
/v reg-8 ^ /v reg-6 



(3.29) 



In Figure [3] we overlay the 3D plots of the rescaled 1Z$ M (m) and lZ^\m) octagon 




Im(m) 



AM (2) 

Figure 3: The superposition of the 3D plots for 1Z S (to) and 1Z 8 (m). 

remainder functions. Specifically, the latter was obtained through the evaluation of 
121 data points using kinematics A, obtained by varying the cross-ratios as 

X + = e<~ 1+2a ' m \ a = 0,...,10, 

X ~ = e n{ - 1+2b/10) , 6 = 0,..., 10. (3.30) 
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Figure 4: A "bird's eye view" of Figure^ computed in kinematics A. 

A similar evaluation was also performed using kinematics C. There appears to be an 
impressive numerical agreement between the rescaled remainder functions at strong 
and weak coupling. 

Figure [4] shows a "bird's eye view" of Figure [3] computed using kinematics A. 
The green regions show where the two-loop rescaled remainder function dominates 
the corresponding strong-coupling result. This could be viewed as an indication of a 
numerical discrepancy between the two. If this is a genuine (rather than purely nu- 
merical) discrepancy, it should be invariant under transformations generated by the 
dihedral symmetry (cyclic permutations and parity) plus conformal transformations]^] 



which manifest themselves 29 as <p — > — (ft and — > <p + ir/2. The green region in 
Figure [4] is clearly invariant under one element of the group, namely the reflection, 
<p — > 7r/2 — <fi. Actually, this symmetry is manifest for our numerical computation in 
kinematics A, one can easily check that it does not involve conformal transformations 
and one is evaluating exactly the same integrals before and after the transformation. 
However, the remaining symmetries do require conformal transformations of the kine- 



7 See Table 1 of 20 , where the regular hexagon case is obtained for all six-point cross ratios equal 
to 1. 

8 Conformal symmetry is necessary because the dihedral group alone transforms the twelve con- 
formal cross-ratios as well as the remaining eight kinematic Lorentz-invariants. To undo the transfor- 
mations on the latter while keeping the former unchanged, one requires conformal transformations. 
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matics and we can see that they do not map the green region to itself exactly. This 
indicates that the green region is not a genuine discrepancy. 




Figure 5: This graph is a plot of 1Z 8 (to) and 1Z 8 (to) as a function of \m\ for 
<fi — tt/4 computed in kinematics C. 




Figure 6: Here we showlZ 8 (m) andlZ 8 (m) as a function of\m\ forcj) = computed 

AM (2) 

in kinematics C. We note that 1Z 8 (m) and 1Z 8 (m) were defined to coincide at 
\m\ = and \m\ = oo. 

We now proceed with a more detailed comparison of our weak-coupling results 
to the strong-coupling remainder function. Specifically, we find it useful to consider 
two-dimensional radial as well as polar plots. In the radial plots, we fix the argument 
of m and vary \m\. In Figure [5] we take a diagonal slice, = 7r/4, whereas in Figure [6] 
we fix the argument of m at = 0. The weak- and strong-coupling results agree 
within the errors - the mismatch between the two is barely visible for both cases. 
By comparing the two radial plots in Figures [5] and [6j we note that the potential 
mismatch between the two depends on the value of <fi. To investigate this further 
we can zoom in on this mismatch at fixed values of \m\, and plot the two remainder 
functions as functions of (ft. 
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Figure |7| shows a polar plot at \m\ 
strong-coupling prediction 7Z 8 (m 



0.2 where we vary between and 7r/2. The 



is depicted alongside our two-loop data 1Z 



(2), 



obtained for the three kinematic sets A, B and C. We notice that the mismatch 



-0.205 



-0.210 



-0.215 



-0.220 



-0.225. 



0.5 



1.0 



1.5 



_ AM (2) 

Figure 7: Angular plot of 1Z S (m) (solid curve) and 1Z 8 (m) for kinematic sets A 
(blue), B (green) and C (red) at \m\ = 0.2 as a function of in the first quadrant. 
The two-loop results include the estimated error-bars. For clarity we have slightly 
shifted the green and red data points. 

— (2) 

between the three conformally equivalent sets for 7Z 8 (m) (which is a good indication 
of the limits of our numerical precision) is comparable with the mismatch between 
the strong- and weak-coupling data, and the remainder functions at strong and weak 
coupling are in agreement within the errors. Figure [8] shows a different polar plot, 
obtained by fixing |m| = 0.4, with similar features compared to Figure [7j 

It is important to note that the overall shift between the curve at strong coupling 
and our perturbative data is extremely sensitive to the precise numerical value for 
the regular octagon TZ^ g _ 8 remainder function. In particular, even the digit in the 
fourth decimal place affects the results in the two polar plots in Figures [5] and |6j For 



these figures, we have used the result quoted in (3.28). 



To conclude this section, we compare the strong- and weak-coupling remainder 
functions for large \m\^ In this region, the integral appearing on the right hand side 
of (3.12) becomes exponentially small, and the strong-coupling remainder function 
approaches the constant value 7tc/6, i.e. twice the result for a hexagon, as discussed 
in Section |3.2 Specifically, we have evaluated the two-loop remainder function at 



\m\ = 1 and \m\ = 1.5, varying the phase from to tt/2 using kinematics A and C. 
The numerical errors in this asymptotic region are larger than the maximal variation 
expected from the strong-coupling side, hence our tests are less stringent here. In 
particular, it becomes difficult to detect the variation in of the two-loop remainder 



3 We would like to thank Fernando Alday and Juan Maldacena for suggesting this additional test. 
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AM (2) 

Figure 8: Angular plot of 1Z 8 (m) (solid curve) and 1Z 8 (m) for kinematic sets A 
(blue), B (green) and C (red) at \m\ = 0.45 as a function of 4>. 



function. The average value we found for in the region is fully consistent with 
twice the result for the two-loop hexagon remainder. 



3.4 Universality of the octagon remainder function 



The agreement between the two shifted and rescaled remainder functions points at a 
linear relation between ( 3.12[ ) and the weak-coupling remainder function, namely 



AM/ 



[U 



13 > 



+ B™ 



(3.3i; 



where Uij denotes the various cross-ratio, and A and B are two constants that we now 



determine. We can evaluate (3.31) at special values of m, in particular for m — 0, 



corresponding to a regular octagon, and for \m\ — > oo, where 7Z$ — » TZ& at strong and 
weak coupling. Doing so, we obtain 



,(2) 
"8-reg 



,(2) 
"6 



tjAM 
/v 8-reg 



2U AM > 



AM-pP) 



8-rce 



-p(2) tjAM 
/v 6 /v 8-reg 



2^ M 



-ip AM 
/v 8-reg 



(3.32) 
(3.33) 



numerical estimates for A^ and B^ 2 \ 
A {2) = -5.41 ±0.03 



Using (|3.19|) and (|3.22|), as well as (|3.27) and (3.28), one arrives at the following 

B {2) = 14.5 ±0.1. (3.34) 



17 



(2) 
reg— 6 - 



Note that the value for quoted above is compatible witfJ 10 -vr 4 /18 = 211 
It would clearly be of great interest if one could find a physicatinterpretation for the 
quantity Af*\ 



An alternative, statistically meaningful way to estimate the values for the param- 
eters and consists in performing a least squares fit of our numerical data 
to the corresponding strong-coupling result calculated using (3.12). We used a set 



of 121 data-points in kinematics A, obtained by varying the cross-ratios according to 



(3.30). From this we obtained the following estimates: 



-5.434 ±0.012 



14.505 ±0.045 



(3.35) 



This can be compared to the numerical values for these parameters in (3.34). The 



values of A^ and in (|3~35|) and KM are in agreement within numerical errors. 
The non-reduced x 2 -square of our fit is r 
of freedom, which in our case is 12 119. 



31, to be compared to the number of degrees 
This is a rather small result, pointing at a 
possible overestimate of our error, still confirming that our data are compatible with 
the strong-coupling result of [29] . In Figure [9] we present the least squares linear fit 
together with our numerical results. 



The errors as quoted by our numerical evaluation routines are usually larger than 
the discrepancy between the different values of the remainder function obtained for a 
fixed set of cross-ratios but with conformally equivalent kinematics. We believe that 
this discrepancy across conformally equivalent sets, which we found to be roughly of 
the order of 10~ 3 , is a slightly more accurate way to quantify (and in fact reduce) our 
errors, although perhaps too optimistic across the whole set of data points. Of course, 
a x 2 -test performed with smaller errors is more stringent. For example, by assuming 
that the numerical error is roughly 1/2 the error as quantified by our routines, the 
new x 2 -square becomes roughly equal to the number of degrees of freedom. 

To summarise, the numerical data we have obtained lead us to conjecture that 
the shape of the remainder function for the octagon Wilson loop is universal. In the 
following section we will further sharpen the statement of shape-invariance of the 
remainder function, and apply it to generic 2n-gon Wilson loops. 

10 In fact, if the errors were actually smaller than what we quote in ( |3.34[ ), the numerical estimate 
A (2) = (-5.5284 + fg)/(-7F log 2 2 + jtj) seems to be in pretty good agreement with -7r 4 /18. 

n We will explain in Section|4]that the linear shift B should actually be absorbed into a redefinition 
of the remainder function at strong coupling. 

12 We recall that the number of degrees of freedom is the number of data points minus number of 
estimated parameters. 



18 



3.670 



3.675 



3.680 



3.685 



Figure 9: The blue line represents the least squares linear fit which maps the strong- 
coupling remainder function (on the horizontal axis) to our weak-coupling results (on 
the vertical axis). The numerical data are also shown in red for comparison. 

4 Remainder functions for general 2n-gons 



Here we will discuss the shape-invariance of the remainder functions in terms of the 
original unrescaled and unshifted remainders for general polygonal Wilson loops (with 
an even number of edges 2n). 

First of all, what do we mean by shape-invariance, or invariance of the functional 
form up to shifts and rescalings? Clearly this invariance (if it persists for general n 
and for general kinematics) implies a linear relation between the remainder function 
at strong coupling and its perturbative expressions at each loop order L: 

^(K'}) = 4n ) < r ° ng ({^}) + , (4.1) 

where and are constants. We assume here that the remainder functions 
at strong and weak coupling have the correct normalisation under simple collinear 



limits, i.e. 18,20 



n 2n TZ 2n -i , (4.2) 



where no constant term can appear on the right hand side of (4.2). By taking simple 



collinear limits of (4.1), one obtains the consistency condition 



A ( L ) strong R (L) _ AL) ^ strong . R (L) ( A o\ 

A 2n /C 2n-1 "+" D 2n ~ /1 2n-l /C 2n-1 + D 2n-\ • K^- 6 ) 
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From this relation, one concludes that A and B must be n-independent. Next, we 



consider triple collinear limits of (4.1), and using (2.12) we arrive at the consistency 
condition 



A ( L ) <r> strong , a (L) ^ 
/1 2n /C 2n-2 + A 2n /C 



strong 



B. 



(i) 

2;i 



j (L) strong 
/1 2n-2 '^2n~2 



<r> >u»n; r(^) _l 7?^) 
/v^„„ 5 -°2n-2 "r '^6 



(4.4) 



Importantly, in the particular two-dimensional kinematics considered in 29 , the final 



six-point configuration has no free cross-ratios left, i.e. the six-point remainder func- 



tion is just a constant. Equation (4.4) implies immediately that A is n-independent, 
in agreement with simple collinear limits. Furthermore, re-expressing TZ^ in terms 



of the corresponding strong-coupling result using (4.1), we arrive at 



B. 



2n 



R (L) , R (L) 
^2n-2 + ^6 



(4.5) 



and hence B 



2n 



[n 



2)B, 



derived earlier that B^ is also n-independent unless B^ 



This result is clearly inconsistent with the result 

0. In turn, this implies 



that B 



2n 



for all n. 



In conclusion, simple and multi-collinear limits show that the weak-coupling result 
for the remainder function and its strong coupling counterpart must be proportional, 



^(K-}) = ^ (i) < rong (K}), 



(4.6) 



where A depends on the coupling constant but not on the number of particles. No 



constant shift is allowed on the right hand side of (4.6) 



We found earlier in (3.31 ) that a nonzero value for B was required in order to match 



the weak-coupling result for the eight-point kinematics and the strong-coupling result 



(3.12) of 29 . There is no contradiction between this and (4.6); indeed, as already 



pointed out in 29 , on the right hand side of (3.12) there is room for an additive 



constant, which was not determined in that paper. Setting 



K 



strong 



n^ M + c 



(4.7) 



and using our earlier result = 
we immediately find C = B^/A^ 



A ( 2 )^strong 



+ 5 (2) , together with (|4~6"]) and ((47} 



Notice also that consistency with multi-collinear limits requires that 

C 



K 



strong 



(4- 



where we have also used that TZ^ M — > 21Z§ M under a triple-collinear limits. 



Finally, we note that, at first sight, it appears that the universality of the re- 
mainder function implies the loss of the maximal transcendentality principle for the 
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(2) 

remainder function TZ 2r [. This is because the strong-coupling octagon remainder func- 
tion (3.12) does not have a uniform degree of transcendentality in the first place, and 
rescalings are not sufficient to restore it. However, one should bear in mind that the 
kinematics we have been investigating is rather special (the remainder function in this 
kinematics lives on a two-dimensional subspace of the full twelve-dimensional space of 
cross-ratios) and this might obscure the transcendentality properties of the remainder 



function. As a concrete example, 13 one may consider the function log(l + x)/a;, which 
appears to have transcendentality equal to zero as x — > 0. This interesting issue 
clearly deserves further study. It would also be very important to investigate the 
validity of the universality conjecture for different number of particles and arbitrary 
kinematics, and confront it with analytic computations at two loops. 



5 Regular 2n-gons 



In this section we will discuss first the particular kinematics associated to regular 
polygons with 2n edges, and we will then move on to summarise the results of our 
numerical analysis we have performed on Wilson loops with up to 30 edges. 



5.1 The regular polygon kinematics 



Here we define the kinematic variables we use for the regular 2n-gon Wilson loop. 
Consider first the projection of the Wilson loop on the (X, Y) plane. This pro- 
jection defines a regular polygon and the positions of the vertices are given by 
Zk = exp(i2irk/2n), where Z := X + iY. In the time direction, the momenta proceed 
following a zig-zag pattern, which leads to the following positions for the vertices (in 
the notation (t, Z))\ 

x 2k = (2sin— , e m — ), x 2k+1 = (0, e"^) . (5.1) 

It is immediate to check that {x 2k — X2k+i) 2 = 0. The expression of the various 
Mandelstam variables needed in the calculation are given by 

2 • 2 ^ 3 

x 2k2j = -4 sin 7T , 

J n 

2 _ 2 

X 2k+12j+l ~ X 2k2j > 

x\ k2j+l = 4(sin 2 |--sin 2 fc ~ J ~ 1/2 7 r) (5.2) 

13 Many thanks to Zvi Bern for pointing out this example and a very useful discussion of how such 
loss of transcendentality for special kinematics occurs in certain two loop integrals. 



21 



At six points, there are only three independent cross-ratios in (2.10) given by 16 19 



«36 



/y>2 ,-v-»2 

x 31 x 46 

rf 2 ™ 2 

x 36 x 41 



f 1 



«14 



/y>2 ,y»2 

x 15 x 24 

/-y» 2 ,-y»2 

x 14 x 25 



"2 



«25 



,->-»2 ™2 
x 26 x 35 
/->-» 2 , 2 
x 25 x 36 



% • (5.3) 



It is immediate to check that for a regular hexagon u\ — u% = 1/3 = 1. For a regular 
octagon, there are four cross-ratios of the type uu+4, and eight of the type uu + z. A 



short calculation shows that, in that case, 



u ii+4: 



1 . 



(5.4) 



Finally, we consider a generic 2n-gon configuration. We find that there are only two 
possible results for the cross-ratios Uij depending on whether % — j is even or odd. 
Our result is 



U 



sin 



sin 



» - J 
i - j 



odd 
2a , 



(5.5) 



where a is an integer. 



5.2 Results at strong and weak coupling 

At strong coupling, Alday and Maldacena in [29] considered the area of the minimal 
surface associated with a regular 2n-gon. In particular they computed the A 



contribution (in the notation of ( 3.14[ )) to the remainder function: 



sinh 
2n 



^ = vr (jn - 2 + I) . (5.6) 

They also considered the n — > 00 limit, where the polygon tends to a circle. In this 
case, the area develops a divergence proportional to the perimeter of the loop, and 
the renormalised area is 

A mn = -2tt . (5.7) 
One may wonder if the corresponding result for the remainder function at weak cou- 



pling shares any of the features of (5.6) 



With the numerical approach outlined in the previous sections, we have computed 
the (complete) remainder function for regular polygons with up to 30 edges at two 
loops in perturbation theory. The regular Wilson loop was calculated in the kinemat- 



ics (5.2), and the remainder function was obtained by subtracting the BDS ansatz. 



The remainder function depends on n through the regular polygon cross-ratios in 
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(5.5). In our perturbative computation we have made no attempt to isolate the A smh 
contribution and our results below will always correspond to the complete remainder 
function. 



Our data are recorded in the Table [2] below, and show quite clearly that, at large 



(2) 

n, 7Z 2 n ~ n reproducing the behaviour of the strong-coupling result (5.6). 



2n 


8 


10 


12 


14 


16 


18 


20 


22 


24 


30 


/c 2n 


-5.528 


-8.386 


-11.262 


-14.145 


-17.035 


-19.926 


-22.821 


-25.717 


-28.614 


-37.311 



Table 2: Numerical values for the Wilson loop remainder function corresponding to a 
2n-gon contour. For a four-edged loop the remainder function vanishes, as there are 
no non-trivial cross-ratios at four points. For the six-point regular polygon, all the 
cross-ratios are equal to one, and the corresponding value of the remainder function 



is given in (3.26), see also (3.27). 



We have also performed a fit of our data to a linear ansatz 1Z 



(2) 

fit,i 



n 



r-x/n. The result of such a fit, using all the data points in Table pi gives 



r x n + r + 



K 



(2) 
fit,l 



7T 



0.9255 n + 2.026 



0.346 



(5i 



Note that n is half the number of edges, In Figure 10 we compare the data points 



with the fit (5.8). The difference between the fit and the data is < 0.01. 




Figure 10: The data points of the regular 2n-gon remainder are well approximated by 
the fit (5.8). The horizontal axis gives the number of edges 2n from 4 to 30. 



In (|5.8|) and (5.9) below we quote the result up to the first digit affected by the error. 
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To see how robust the fit is, we have also considered a slightly more general fit 
with an additional 1/n 2 term, and found 

( 2 ) ( 0.181 0.228\ , , 

^fit,2 = vr(^- 0.9239n + 1.9955 —J . (5.9) 

The difference between this fit and the data is < 0.0004. One could also perform 
fits including higher inverse powers of n, however these do not change our conclusion 
about the linear behaviour of the data as a function of n at large n and the estimate 
for the slope of the line. 

A few comments are in order. Firstly, we note that at large n the n-dependence 
of is the same as that of A smh at strong coupling - they both grow linearly in 
n. Secondly, we observe that the constant term is well approximated by 2ti and up 



to a sign is the same as the strong coupling result in (5.6). We have already noticed 



earlier a similar sign swap between strong and weak coupling results, see Figure [2j 

Here we will not attempt to make a precise connection between these results 
and the universality of the functional form of the remainder function advocated in 
earlier sections. To do so, we would require the full expression for the remainder 



function at strong coupling. In addition to A smh , given in (5.6), we would also need 
expressions for AAf° s and A^ Ia (which can be obtained in principle from [29] but 
are not known explicitly for arbitrary n). In addition to these, as already noted 
in 29 , there is a regularisation-prescription dependent shift. This shift is a kinematic- 
independent constant, which does depend on n. Therefore, for regular polygons it 
is not immediately straightforward to give a direct relation between the complete 
remainder functions at strong and weak coupling. 



6 Conclusions 



In this paper we have presented strong evidence in favour of a remarkable albeit 
unexpected agreement between the remainder functions computed at strong and at 
weak coupling. We would like to comment briefly on the conclusions one can draw 
from this. 

1. Our evaluation of the two-loop remainder function is mostly numerical. There- 
fore there is always the possibility that there exists a genuine disagreement between 
the strong- and weak-coupling remainder functions. If this is the case, it is still 
remarkable that it must be very small, as quantified in this paper. 

2. The second possibility is of course that there is a precise agreement between 
the remainder functions at strong and weak coupling. Even if this holds only for eight 
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points, and for the special kinematics of 29 , this would be a highly non-trivial result 



but it is more natural to conjecture that it holds generally for any number of points, 
and for arbitrary kinematics. In this case, the remainder function could be cast in 
the form 

K(a) = A(a)K stmng , (6.1) 

where the function A (a) has a perturbative expansion starting at two loops and is 
independent of n and of the kinematic variables (which are all contained in 7^ strong ). 



If (6.1 ) is correct, it would be of great importance to understand the physical meaning 



of the quantity A (a) 



The conjecture of universality for the Wilson loop in (6.1 ), would imply novel non- 
renormalisation theorems or symmetry arguments. An exciting way forward would be 
to develop analytic methods for computing Wilson loops at strong coupling for other 



values of n and/or in different kinematics, perhaps by extending the results of 29 
Our conjecture of universality of the Wilson loop remainder function could then be 
further tested in these cases. 
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